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DESIGN  PROCEDURE  FQt  TRANSPIRATION  COOIED 
GAS-TURBINE  BIADES 


Summary 


Contained  herein  is  the  calculation  procedure  for  determining  the 
coolant  flow  required  to  maintain  a  given  wall  temperature  in  transpiration* 
cooled  turbine  blades.  The  basic  method  is  not  new;  it  has  evolved  from 
numerous  modifications  based  upon  the  experience  acquired  during  testing  and 
evaluation  of  transpiration  cooled  configurations.  The  theory  supporting 
the  equations  is  included. 

The  calculation  procedure  presented  in  this  report  has  been  tried 
and  has  proven  to  be  a  useful  tool.  It  is  believed  that  continued  testing 
will  further  the  understanding  of  the  influence  of  turbulent  flow  on  the 
aerodynamic  and  heat  transfer  characteristics  of  this  type  of  cooling.  With 
refinement  of  the  theoretical  approach,  improvement  in  the  design  procedure 
can  be  realized. 

In  order  that  potentially  feasible  configurations  can  be  evaluated 
quickly  and,  therefore,  economically,  the  calculation  method  has  been  pro¬ 
grammed  for  the  IBM  70u  digital  computor.  A  brief  description  of  the  methods 
of  preparing  input  data  for  the  machine  and  of  the  operation  sequence  employed 
by  the  machine  is  also  included. 
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Object 


1.  To  outline  in  detail  the  procedure  used  in  the  design  of  transpiration 
cooled  turbine  blades. 

2.  To  evaluate  the  results  of  cascade  testing  for  the  purpose  of  confirming 
the  theory  supporting  the  equations  in  the  calculation  procedure  and  to 
incorporate  those  modifications  deemed  necessary  for  improving  the  accuracy 
e  *  the  calculation  method, 

3.  To  program  the  refined  design  procedure  into  the  I.B.M.  7 Ok  Computer. 


Conclusions 


1.  A  refined  design  procedure  has  been  worked  out  in  detail  to  aid  in  the 
design  of  transpiration  cooled  turbine  blades. 

2.  The  refined  design  procedure  has  been  programmed  on  the  70U  Computer, 

thus  enabling  the  evaluation  of  many  cooling  configurations  heretofore 
considered  impractical  because  of  time  and  cost  limitations. 


Re: ommendations 


10  Employ  the  procedures  outlined  herein  in  the  design  of  transpiration  air 
cooled  turbine  blades, 

2.  Continue  cascade  testing  to  farther  the  correlation  between  theory  and 
experimental  data. 

3.  Utilize  the  I.B.M.  70U  Computer  program  developed  herein  to  facilitate 
the  design  of  transpiration  cooled  configurations. 
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SYMBOLS 

The  following  symbols  are  used  In  this  reports 
A  Cross-sectional  area,  sq.  ft. 

a.  Coefficient  in  linear  spanwise  variation  of  cooling -air  temperature,  Tc*Tc^r  «3L x 

5  Chordwise  peripheral  width  of  coolant  passage,  ft 

b*  Skin  area  associated  with  a  particular  passage  and  cross  section, 

/3  Angle  between  blade  axis  and  rotor  radius  sector,  degrees 
Cp  Specific  heat  at  constant  pressure, 

£  Leading  edge  diameter,  ftj  D*  Hydraulic  diameter- of  coolant  passage,  ft 
Deiivitive 

F, .  Euler  number  of  gas  flow  over  blade, 

w  o  y 

f  Friction  force,  lbs 

f  Coolant  passage  friction  fac tori  vw  Coolant  flew  parameter  for  laminar  flow  region 

0  Acceleration  due  to  gravity,  -F+  /see  * 

h  Heat-transfer  coefficient,  - — a  o 

hr*  v-t  K 

K  Permeability  coefficient,  sq.  ft. 

k  Thermal  conductivity,  - 

•  v  I'  'SttC 

L  Length  from  stagnation  point  to  trailing  edge  along  either  surface,  ft 
m  Mass  coolant  flow,  Lb  /s<?C 


My  Nusselt  number,  h  I  /  k 

_  ^  C  p 

CO  Angular  velocity,  sec  “j  Xl r — ^ — 

p  Pressure 
Pr  Prandtl  number,  c 

K 

Density ,,  Lb /-ft* 

^  Heat  transfer  by  convection  tc  blade  skin  from  hot 'gas 

/ 
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q^v  Hsat  transfer  by  radiation 

c^s  Heat  transfer  by  sidewise  convection  in  the  coolant  boundary  layer 
^  Heat  transfer  by  conduction  through  blade  wall 
Re  Reynolds  number  defined  by 

/^w 

Reynolds  number  defined  by 


w 

r  Radius,  ft,  r*  « 

R  Gas  constant,  -Tt  -  tb/tb°  F 


T  Temperature  °R 

T  Porous  blade  -  shell  thickness,  ft 
Tad  Adiabatic  gas  temperature,  °R 

U  Free  stream  gas  approach  velocity  upstream  of  blade  leading  edge,  ft  /sec 
‘J  *  Dimensionless  roans  velocity  in  fl*ea  stream. 

3  *  72*  ft  _ 

(  o 

\r  Cooling  air  velocity  through  porous  wall,  ft /see 

/.*-  Absolute  viscosity  (based  on  porous  wall  temperature).  Lb /-ft.  50c 

W  Gas  velocity  relative  to  blade,  -TL/scc 

Wi  Coolant  velocity  in  passare,  ft  /  sac, 

X  Spanwise  distance  in  coolant  passage  from  blade  root,  ft 
y  Peripheral  distance  from  stagnation  point  along  either  surface  of  blade,  ft 
y  *  Dimensionless  distance  from  stagnation  point  along  blade,  y/L 
1  Distance  normal  to  skin  in  the  gas  direction,  ft 
Subscripts: 

C  Cooling  air  or  cooling  air  passage 
Denotes  coolant  into  air  passage 
C mC.  Denotes  coolant  out  of  air  passage  1 
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Z  External  gas  flow 
4  Internal  coolant  flow 

q  Refers  to  fixed  point  in  stream  (e.g,,*  upstream  condition) 
r  Coolant  passage  entrance  (blade  root) 
jj  Stream 
w  Porous  wall 
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Introduction 


As  analytical  theory  has  predicted,,  testing  has  shewn  that  trans¬ 
piration  cooling  is  a  very  effective  method  of  cooling  gas  turbine  blades. 

This  report-  intends  to  provide  a  better  understanding  of  the  transpiration 
cooling  process  by  summarizing  the  design  procedure  that  has  evolved  during 
five  years  of  usage. 

Inasmuch  as  porous  materials  currently  available  for  use  in  trans¬ 
piration  cooled  turbine  blades  do  not  possess  sufficient  strength  to  with¬ 
stand  the  high  stresses  imposed  by  rotative  turbine  speeds,  an  internal  load 
carrying  ’•  strut”  is  required  to  which  the  porous  skin  is  attached.  In  ad¬ 
dition  to  serving  as  the  blade  supporting  member,  the  strut  divides  the  cross- 
section  of  the  blade  into  compartments. 

Because  of  the  acceleration  of  the  gasses  through  the  rotor,  there 
exists  around  the  periphery  of  an  individual  blade  a  substantial  variation  of 
static  pressure.  The  amount  of  coolant  flow  that  passes  through  the  porous 
wall  is  dependent  on  the  difference  in  the  squares  of  the  absolute  static 
pressure  levels  on  opposite  sides  of  the  wall  with  the  result  that  small  changes 
in  pressure  level  can  cause  very  large  changes  in  coolant  flow  rate.  It  has 
also  been  observed  that  the  temperature  distribution  changes  drastically  as  the 
coolant  flow  rate  is  varied.  Thus,  the  difficulty  in  achieving  uniform  temper¬ 
ature  around  the  periphery  of  a  blade  is  readily  seen. 

The  optimum  design  for  a  transpiration-cooled,  strut -supported  turbine 
rotor  blade  as  regards  coolant  flow,  therefore,  incorporates  a  blade  skin  with 
a  variable  chordwise  and  a  variable  spamd.se  permeability.  This  is  attributed 
mainly  to  the  fact  that  the  large  variations  in  external  gas  pressure  and 
temperature  give  rise  to  similarly  large  variations  in  coolant  demand.  At  a 
specific  chordfod.se  position,  for  example,  the  midsection  of  the  blade  may  "see" 
a  gas  temperature  200°F  higher  than  the  tip  section,  and  to  make  matters  worse, 
a  higher  static  pressure  level.  Thus,  for  a  given  coolant  pressure  level,  a 
material  with  constant  permeability  would  have  to  be  chosen  to  satisfy  the 
maximum  coolant  requirement  with  consequent  wastes  of  coolant  in  the  areas  of 
lesser  demand.  As  the  peripheral  size  of  the  passages  decreases,  these  effects 
in  the  chordwise  direction  become  less  poignant  except  in  the  leading  and 
trailing  edge  areas. 

The  problems  incurred  in  the  fabrication  of  blade  skins  with  variable 
permeability  are  still  found  prohibitive.  Consequently,  the  underlying  princi¬ 
ple  in  the  procedure  is  to  assume  a  permeability  and  coolant  supply  pressure 
level  which  will  yield  satisfactory  wall  temperatures. 


The  design  method  presented  herein  prescribes  a  uniform  skin  temper¬ 
ature  equal  to  the  maximum  allowable  for  the  material.  In  this  way,  the  coolant 
demand  can  be  held  to  a  minimum.  The  division  of  the  blade  interior  into 
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separate  passages  and  the  use  of  sized  orifices  in  the  blade  base  to  meter 
coolant  to  the  various  passages  tend  to  minimize  local  overcooling  and  total 
bleed  requirement. 

This  report  attempts  to  set  down  an  outline  of  the  basic  principles 
of  transpiration  air-cooled  blade  design.  It  deals  only  with  the  thermodynamic 
aspect  of  cooling  methods.  No  consideration  is  given  to  the  important  problems 
associated  with  stress  analysis  and  fabrication.  Moreover,  the  effects  of 
changes  in  engine  speed,  flight  speed,  and  altitude  on  pressure  levels  and 
pressure  distributions  around  the  periphery  of  turbine  blades  is  not  discussed. 
However,  it  should  be  stated  here  that  inasmuch  as  the  permeability  of  trans¬ 
piration  cooled  blades  must  be  designed  to  provide  the  proper  cooling  at  a  pre¬ 
scribed  pressure  distribution,  variations  of  flight  speed,  engine  speed,  and 
altitude  may  cause  the  permeability  to  be  inconsistent  with  the  aerodynamic 
pressure  distribution.  The  uniform  design  temperature  will  no  longer  prevail 
and  in  the  process  of  providing  adequate  coolant  to  the  hottest  portion  of  the 
blade,  excess  coolant  flow  will  result  in  the  cooler  areas; 


t 
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Results 


As  progress  is  made  in  predicting  convection  heat  transfer  coeffi¬ 
cients,  the  design  procedure  outlined  in  this  report  can  be  steadily  improved* 
The  reliability  of  any  such  method,  however,  can  be  established  only  through 
comparison  of  the  results  of  well  instrumented  tests  with  predicted  values* 

In  this  way,  modifications  of  the  theory  can  be  made  so  that  a  fairly  good 
design  method  is  assured* 

Hie  cascade  data  developed  under  Contract  NOas  5>6-b9!>**c  (Reference  8) 
should  be  compared  with  the  results  predicted  by  the  method  presented  herein* 
It  is  unlikely  that  agreement  will  be  exact  in  light  of  the  many  assumptions 
made  in  the  analytical  approach  and  the  errors  inherent  in  the  experimental 
work,  but  a  study  of  this  kind  can  aid  in  establishing  the  limitations  of  the 
present  approach* 

Since  computations  can  be  made  rapidly  on  the  70b  computor,  some 
interesting  design  studies  can  be  made  which  could  hardly  be  justified  if 
manual  calculation  were  employed.  For  instance,  the  chordwise  temperature 
variation  along  the  skin  of  coolant  passages  in  regions  where  the  external  gas 
pressure  changes  rapidly  is  of  interest.  The  most  important  example  is  the 
leading  edge  passage,  Ihe  calculations  can  be  made  by  dividing  the  passage 
into  a  number  of  fictitious  passages  with  identical  coolant  pressure  distribu¬ 
tions*  Actually,  this  is  the  only  way  to  determine  the  total  coolant  require¬ 
ment  for  such  a  passage. 

The  analytical  determination  of  the  heat  that  must  be  removed  from 
the  blade  skin  at  the  stagnation  point  along  the  leading  edge  does  not  take 
into  consideration  the  heat  conduction  within  the  skin  to  cooler  areas  (on 
the  pressure  and  suction  sides).  This  safety  factor  will  have  an  equalizing 
effect  on  the  chordwise  distribution  of  the  blade  skin  temperature, 

A  resume  of  the  design  procedure  as  described  herein  follows! 

1,  Assume  a  value  for  the  internal  static  pressure  at  the  root  section 
(P5,r)  and  calculate  the  static  pressure  at  any  spanwise  position 
(Ei,x)  from  equation  (5)f 

2,  Assume  a  value  for  the  blade  wall  temperature  (Tw)  and  calculate 
inside  gas  velocity  (v)  from  equation  (9), 

3,  Calculate  coolant  flow  parameter  (fw)  from  equation  (19)  using 
the  calculated  value  of  (v),  (A  modified  equation  for  (f^)  i« 
used  at  the  leading  edge  (l)). 
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b.  Calculate  Euler  number  (Eu)  and  the  ratio  of  adiabatic  gas  tempera 
ture  to  wall  temperature  (T^/  • 

5*  For  the  laminar  flow  region,  obtain  -j  '  from  Figure  (h),  and 
calculate  Tw .  For  the  turbulent  regfbn  use  equation  (2i*)  to  calcu¬ 
late  . 

O,  If  the  calculated  value  of  Tw  does  not  check  the  assumed  value, 
the  calculation  is  repeated  for  a  new  •  If  the  values  obtained 
for  Tw  are  now  compatable  with  design  requirements,  a  new  p;  r 
should  be  tried.  ’ 

7.  "When  a  satisfactory  skin  temperature  has  been  obtained  at  all  loca¬ 
tions  for  which  a  calculation  has  been  made,  the  mass  flow  (  /&\r 
Lb./s^C.  Ft2)  is  determined  at  each  location  from  equation  (9). 

8.  The  coolant  requirement  for  each  passage,  rrv  is  obtained  by  sum¬ 
ming  the  product  of  each  /°V*  and  associated  skin  area  over  the 
entire  passage.  The  sum  of  the  mg  quantities  for  all  passages 
represents  the  total  coolant  requirement  per  blade. 

9.  The  equations  required  to  obtain  the  orifice  area  necessary  for  tht 

specified  Pj  are  contained  in  reference  (l). 
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Discussion 


Statement  of  the  Design  Problem 

For  a  specified  set  of  turbine  design  parameters,  i.e.  rotor  speed, 
gas  temperature  and  pressure,  blade  geometry,  coolant  temperature  and  pressure, 
etc.,  the  temperature  of  the  transpiration  cooled  blade  skin  will  depend  upon 
the  porosity  distribution  of  the  skin  and  the  pressure  distribution  in  the 
coolant  passages.  The  problem  can  be  formulated  in  three  ways: 

I.  Determine  the  required  coolant  pressure  distribution  to  give  a  uni¬ 
form  wall  temperature  with  a  specified  porosity  distribution. 

II.  Determine  the  porosity  distribution  required  with  a  specified  coolant 
pressure  distribution  to  give  a  constant  wall  temperature. 

III.  Determine  the  wall  temperature  distribution  resulting  from  specified 
pressure  and  porosity  distributions. 

These  three  cases  are  summarized  in  Table  I* 

Design  Aspects 

Inasmuch  as  the  coolant  pressure  distribution  cannot  exceed  the  com¬ 
pressor  outlet  pressure,  and  must  increase  from  root  to  tip  in  a  rotor  blade, 
it  cannot  be  varied  indiscriminately;  hence.  Case  I  is  not  a  very  practical 
way  to  proceed.  Manufacturers  of  porous  material  have  yet  to  provide  a  con¬ 
vincing  demonstration  of  their  ability  to  provide  a  well  controlled  variation 
in  porosity.  Accordingly,  Case  II  is  not  yet  practical.  The  best  approach 
appears  to  be  defined  by  Case  III.  The  porosity  may  be  specified  constant 
and  the  temperature  distribution  determined  for  a  particular  pressure  distri¬ 
bution  which  is  possible  to  obtain.  If  the  temperature  so  calculated  is  un¬ 
satisfactory,  a  new  pressure  distribution  can  be  determined  and  the  process 
repeated. 


| 

Jt 
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If  a  satisfactory  skin  temperature  cannot  be  obtained  with  the  maxi¬ 
mum  coolant  pressure  available,  at  some  available  porosity,  consideration  may 
be  given  to  adding  an  auxiliary  compressor  to  increase  the  coolant  pressure 
and/or  to  including  a  ram-air  cooler  to  reduce  the  coolant  temperature. 

Analysis 


The  pressure  distribution  in  the  blade  passages  is  calculated  after 
assuming  a  value  for  the  coolant  pressure  at  the  blade  root.  The  root  pressure 
can  be  made  less  than  the  coolant  supply  pres  sure -by  the  insertion  of  an  ori¬ 
fice  in  the  blade  root. 
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The  pressure  along  the  passages  of  stator  blades  will  remain  un¬ 
changed  even  at  relatively  high  flow  rates.  This  was  verified  by  the  cascade 
test  of  stator  blades  developed  under  NOas  $6-k9$-c  (Ref.  8)*  Rotor  blades 
differ,  however,  since  the  centrifugal  force  must  be  considered.  The  calcu¬ 
lation  method  according  to  Reference  3  is  reviewed  briefly  here. 

Consider  a  rotor  blade  channel  of  varying  cross  sectional  area  as 
illustrated  in  Figure  1. 

Assume  steady  flow  through  the  channel,  and  write  Newton’s  second 
law  for  the  fluid  in  the  cross  hatched  volume  element  of  the  channel.  This 
law  must  be  applied  to  a  control  volume  that  is  at  rest  relative  to  the 
channel. 

Newton’s  second  law  states  that  the  summation  of  forces  acting  on  a 
body  in  a  given  direction  must  equal  the  rate  of  change  in  momentum  of  the 
body  in  this  direction.  In  this  case  it  is  applied  to  the  flow  direction. 

The  sura  of  all  pressure  forces  is: 


Pi  A**Cp;  pi  )(A-dA)~  (p,-  +f*Pi)dA 


The  third  term  comes  about  because  the  pressure  at  the  wall  inclined 
to  the  flow  direction  has  a  component  in  the  flow  direction  which  is  approxi¬ 
mately  equivalent  to  the  average  pressure  on  the  wall  (p|’  )  multiplied 

by  the  wall  area  and  the  sine  of  the  angle  ,5,  The  wall  area  multiplied  by  sine 
(3  is  equal  to  <JA.  Neglecting  second  order  terms,  the  total  pressure  force 
becomes  “A  d  Pi  • 

The  friction  force  can  be  expressed  by  the  Fanning  Equation  as 

,  -frWi2Adx 
F  ' - ^ - 

The  body  force,  caused  by  centrifugal  acceleration  of  the  fluid 
particles,  can  be  written  as 

Ar  CO2  Cos  (3  dx 

The  change  in  momentum  of  the  stream  flowing  through  the  control 

volume  is 

drn  Wi  s  W i  +  Wi* dm 


yields 


Equating  the  summation  of  the  force  terms  to  the  change  in  morsertum 


-Adp,f^,ArtJ2Cos.(3dx-^^id^  -  mdWt+Widm  (1) 
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The  number  of  variables  in  Equation  (1)  can  be  reduced  by  use  of  the 
continuity  equation 


m  *  (>  AWI 

to  eliminate  W^ ,  and  the  equation  of  state 

h*  zp  RT 


(2) 


(3) 


to  eliminate  o  .  Additional  simplification  results  from  introducing  the 
relation  /m, 

dm  Ad  CXJ  ,  m  dA 

4i - 4*  X  dx  (W 

As  a  result  of  these  substitutions,  Equation  (l)  becomes 


($) 


In  Equation  (5),  the  specific  mass  flow  m/A  ,  the  Internal  pressure 
Pi,-  the  radius  r,  the  cross  sectional  area  A,  and  the  temperature  T  are  gener¬ 
ally  functions  of  the  distance  x  from  the  channel  entrance  measured  along  the 
channel  axis.  The  x  dependence  of  A  and  r  is  easily  determined.  The  coolant 
temperature  is  assumed  to  increase  linearly  according  to 

T*  Tr-»ax  (6) 

Where  11  an  is  an  estimated  quantity.  Flow  velocities  are  usually  low  enough 
in  coolant  passages  so  that  the  total  temperature  can  be  used  instead  of  the 
static  temperature. 

A  second  expression  including  the  same  variables  as  functions  of 
the  distance  from  the  passage  entrance  can  be  obtained  from  a  consideration 
of  the  pressure  drop  through  the  porous  wall  when  the  pressure  distribution 
along  the  outside  of  the  passage  is  prescribed.  The  mass  velocity  of  the 
coolant  ejected  through  the  porous  wall  at  any  location  multiplied  by  the 
width  of  the  porous  portion  of  the  channel  wall  measured  in  the  circumferential 
direction  equals  the  rate  of  change  of  the  coolant  mass  flow  passing  through 
the  passage  at  that  location;  that  is, 

b  Ipsr)  =  — stra. 

^  '  d*  (7) 

The  quantity  py  is  the  average  mass-flow  rate  over  the  channel  width  b.  The 
negative  sign  appears  in  the  right  member  of  equation  (7)  since  m  decreases 
with  increasing  distance  along  the  coolant  passage.  -  The  mass  velocity  of  a 
gas  flowing  through  a  porous  skin  is  related  to  the  pressure  drop  by 
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Pi*-P«*  .  C,  <2  RT<**)  (e sr)  f  c2  2  RT  ({V)* 
T 


No  general  relationship  for  the  parameters  and  C2  has  been  obtained  as  yet 
so  that  thev  have  to  be  determined  by  experiments  for  each  porous  material. 
Equation  (8)  can  be  simplified  to  / 


<-v*cK(pi2-pe2)r' 


(9) 


where 

Ck--3.0SO(  ZJtJ*  (10) 

For  a  typical  wire  cloth,  n  *  5/8.  Substituting  equations  (1*)  and  (9)  into 
equation  7  gives  the  desired  second  relation: 

A  d  (irhjg  dA  =  -bCt,(P|*-bexr  (11) 

Equation  (5)  and  (11)  must  be  solved  simultaneously  to  yield  p^  (x). 
When  the  local  coolant  rate  is  prescribed,  a  direct  solution  is  possible. 

However,  the  more  important  case  in  which  the  local  coolant  rate  is  a  quantity 
to  be  determined  requires  a  trial  and  error  solution. 

The  iteration  process  can  be  avoided  by  an  approximate  approach 
which  assumes  the  friction  force  and  the  rate  of  change  of  momentum  are  negli¬ 
gible,  That  such  an  assumption  is  justified  can  perhaps  be  better  understood 
through  consideration  of  an  analogy  to  the  present  case. 

If  a  vertical  standpipe  is  filled  with  a  liquid,  a  significant  pressure 
increase  is  present  from  top  to  bottom  because  of  the  body  force  exerted  by  the 
dense  liquid.  If  a  standpipe  were  filled  with  gas,  the  pressure  variation  would 
be  negligible.  However,  if  the  pipe  were  rotated  at  turbine  speeds,  the  body 
force  on  the  gas  will  increase  to  such  an  extent  that  a  pressure  gradient  simi¬ 
lar  to  that  in  a  liquid  exists. 

Now  suppose  that  holes  are  punched  in  the  standpipe  containing  the 
liquid  and  that  a  supply  of  liquid  is  available  at  the  top  of  the  pipe  suf¬ 
ficient  to  make  up  for  liquid  lost  through  the  holes.  If  the  flow  is  not 
large,  the  pressure  distribution  will  not  be  much  different  than  it  was  with 
the  holes  absent.  As  the  flow  rate  is  increased  either  by  adding  more  holes 
and/or  increasing  the  supply  pressure,  the  velocity  of  the  fluid  in  the  pipe 
becomes  large  so  that  the  friction  and  momentum  terms  in  equation  (1)  must  be 
included. 
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The  simplified  equation  becomes  a  force  balance, 

“A  dpj  +  Ara?’dxaO  (12) 

which  can  be  directly  integrated  to  give: 

I "  (-*§?)  =  TEa*  ^'0  +  (M) 

Calculations  using  both  the  approximate  and  the  rigorous  approach  for  rotor 
blade  cooling  passages  show  that  large  errors  are  not  introduced  by  the  ap¬ 
proximate  method  even  at  large  coolant  flow  rates. 

The  coolant  mass  flow  rate  required  to  give  a  specified  local  value 
of  the  skin  temperature  cannot  be  obtained  directly  from  equation  (9)  since 
the  skin  temperature,  Tw,  which  appears  in  this  equation  depends  upon  the 
relative  rates  at  which  heat  is  added  by  the  gas  and  removed  by  the  coolant. 
Therefore,  another  equation  relating  Tw  and  /^v  must  be  obtained  from  a  heat 
balance  at  the  position  in  question. 

A  section  of  the  porous  skin  is  illustrated  in  Figure  2.  Heat  can 
enter  the  bounded  region  by  convection  from  the  hot  gas,  q  :  from  the  entering 
coolant,  i  i  hy  conduction  along  the  porous  wall,  qw;  from  sidewise  con¬ 
vection  in  the  coolant  boundary  layer,  qs;  and  by  radiation,  q^  Heat  leaves 
the  bounded  section  in  the  exit  coolant,  qc.  Heat  transfer  to  turbine  blades 
by  radiation  is  usually  negligible,  and  if  the  sidewise  heat  transfer  terms 
are  also  neglected,  a  heat  balance  on  the  bounded  section  becomes: 


*■>  (Tad  *  Tw)  --Cp(«wVr(Tw-Tfc)dA 


(U») 


The  first  term  represents  the  convection  heat  transfer  from  the  gas  and  the 
term  on  the  right  the  net  heat  picked  up  by  the  coolant.  Equation  (lU)  can 
be  rewritten  in  the  forms 


/O  ^r*  hC"kd  -Tw) 


(15) 


In  order  to  use  equation  (1£>)  in  conjunction  with  equation  (9)  to 
determine  the  coolant  flow  and  Tw,  and  expression  for  the  heat  transfer  coef¬ 
ficient  h  is  the  most  difficult  problem  of  the  entire  analysis. 

The  heat  transfer  coefficient  depends  upon  the  character  of  the 
boundary  layer  and  therefore  will  change  markedly  when  the  boundary  layer 
changes,  from  laminar  to  turbulent.  The  heat  flow  from  the  gas  boundary  layer 
on  the  outside  of  the  blade  into  the  blade  surface  is  given  by 
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k  (^y)w  dA  *  h  (Tad  -  Tw)d A 


(16) 


Equation  (16)  can  be  solved  for  h  to  give:  / 


(17) 


which  3hows  h  to  be  a  function  of  the  temperature  gradient  at  the  vrall  and 
the  difference  between  the  gas  and  wall  temperatures. 


Heat  Transfer  Through  Laminar  Boundary  Layer 


For  the  laminar  region,  an  exact  solution  of  the  boundary  layer 
equation  can  be  obtained  io  yield  •  This  exact  approach  is  difficult 

and  tedious.  The  uncertainties  which  till  are  connected  with  the  calculation 
of  heat  transfer  on  turbine  blades  at  present  do  not  justify  the  time  expendi¬ 
ture  necessary  for  exact  calculations  and  suggest  approaches  which  derive  ap¬ 
proximate  heat-transfer  parameters  with  a  small  effort. 


In  Reference  f>  the  calculation  of  local  heat -transfer  coefficients 
from  approximate  solutions  to  the  laminar  boundary  layer  equations  around  a 
porous  wedge  was  shown.  From  these  results  for  porous  wedges,  local  heat- 
transfer  coefficients  for  the  laminar  region  of  a  gas  turbine  blade  can  be 
obtained  by  the  stipulation  that  the  local  heat-transfer  coefficient  at  a 
specified  point  along  the  blade  periphery  is  identical  to  that  on  a  wedge  for 
which,,  at  the  same  distance  from  the  stagnation  point,  the  Euler  number  is 
the  same  as  that  on  the  blade.  The  Euler  number  is  given  by: 


r  .  r.  %  w 

■a  -  u  a  <  08) 

The  i-esults  of  Reference  £  calculations  are  plotted  in  Figure  3 
which  presents  a  relation  between  the  heat-transfer  coefficient  contained  in 
Iluj  the  coolant  velocity  contained  in  fw,  and  the  known  flow  parameters. 

In  this  figure  the  external  and  coolant  flows  are  assumed  to  have  equal  fluid 
properties  and  $  calculated  for  small  mach  numbers.  The  coolant 

flow  parameter,  is  given  by: 


(u>) 


The  magnitude  of  -fw  is  an  indication  of  the  rate  of -coolant  flow. 

f 
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In  the  discussion  thus  far,  the  following  relationships  have  been 
established; 

(a.)  pj  =  f  (x)  Prom  (5) 

(b) 

Prom  (9) 

(c.)  Pj  *  ^(h,Tw)  From  (15) 

(d)  b  *  From  (Fig.  13) 

This  system  of  equations  contains  four  unknowns  pi,  Tw,  /  ,  h;  and  can  be 
solved  for  the  desired  quantities,  Tw  and  v  ,  Since  the  relation  in  (d)  is 
presented  in  graphical  form,  the  solution  is  a  trial  and  error  process* 

In  order  to  facilitate  calculations.  Figure  U  was  prepared  which 
eliminates  h  from  equations  (c)  and  (d)0  The  calculation  procedure  requires 
evaluation  of  p^  from  (a)  and  an  iteration  process  between  (b)  and  Figure  U* 

Heat  Transfer  Through  Turbulent  Boundary  Layers 

No  exact  solutions  exist  for  the  case  in  which  the  boundary  layer  is 
turbulent.  Experimental  information  is  limited  on  the  effect  of  pressure  vari¬ 
ations.  It  is  known,  however,  that  pressure  variations  have  less  influence  in 
the  turbulent  flow  region  on  solid  surfaces  than  in  the  laminar  flow  region* 

For  flow  with  constant  pressure,  an  approximate  method  developed  in  Reference 
6  and  later  simplified  by  Reference  7  checks  experimental  results  fairly  well* 
The  assumptions  made  in  Reference  7  are: 

le  Tne  flow  consists  of  a  turbulent  region  and  a  laminar  sublayer  which 
separates  the  turbulent  flow  from  the  wall  surface* 

2.  The  temperature  and  velocities  in  the  turbulent  region  have  the  same 
values  on  a  transpiration  cooled  wall  as  in  an  ordinary  boundary 
layer  on  a  solid  surface  under  otherwise  identical  conditions* 

By  restricting  the  investigation  to  fluids  with  a  Prandtl  number  near 
unity,  the  following  expression  for  the  wall  temperature  is  obtained  in  Refer¬ 
ence  7? 

Tw -Tc  r _  rl _ 

2  (20) 

where  r*  is  the  ratio  of  the  velocity  parallel  to  the  surface  at  the  border 
between  the  laminar  sublayer  and  the  turbulent  part  of  the  boundary  layer,  to 
the  stream  velocity  outside  the  boundary  layer  and  where 

si » jihcsSe. 

h 


(21) 


(al  Euler  Number,  1,C'.  (b)  Euler  Number,  0.5«  (c)  Euler  Number 
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with  h  the  heat-transfer  coefficient  that  would  be  present  on  a  solid  surface 
under  the  sans  outside  flow  conditions.  The  equation  for  turbulent  heat  trans¬ 
fer  on  a  solid  flat  plate  can  be  used  in  determining  h. 

* 

Under  these  conditions 
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Moreover.  for  a  flat  plate 


r'«  * 


.JLil 


c<?^  (tj  U'-otw) 


(22) 


(23) 


Substitution  of  equations  (22)  and  (23)  into  equation  (20)  and  simplification 
yields 

2  »f 


7. 

T 


»?.<  p» 


(2I4) 


Equations  (a)  and  (b)  which  were  used  to  obtain  a  solution  for  the 
laminar  flow  region  are  also  used  for  the  region  over  which  the  boundary  layer 
-is  turbulent.  However,  equations  (z)  and  (d).  for  which  Figure  (U)  was  substi¬ 
tuted  in  the  laminar  flow  case-,  are  now  replaced  by  equation  2U. 


The  Com-put or  Program 

Because  of  the  tedious  and  time  consuming  iterative  process  involved 
in  the  calculations  through  equation  (7 )P  this  portion  of  the  method  has  been 
programmed  for  the  ?0U  computor.  In  order  to  establish  a  systematic  scheme  for 
the  machine  calculations,  the  coolant  requirements  and  skin  temperatures  for 
all  of  the  passages  in  each  blade  section  are  evaluated  before  going  to  the 
next  section.  At  each  section  the  calculations  are  made  for  the  leading  edge, 
the  laminar  region,,  and  the  turbulent  region  respectively* 

This  problem  was  prepared  for  the  IBM  70U  digital  computor  using  the 
FORTRAN  system. 

For  a  given  set  of  inputs  corresponding  to  a  particular  airfoil  sec¬ 
tion,,  the  computor  program  begins  with  the  leading  edge  calculations*  Utiliza¬ 
tion  of  a  bivarient  map  is  made  which  given  (1  -  Q>  )  as  a  function  of  fw  and 
T<,/Tv,  assuming  Eu  «  1.0.  This  map  was  tabulated  so  that  the  program  could 
obtain  values  of  (l  «  Cp)  by  interpolation.  The  order  of  interpolation  on  *w 
was  2  and  on  Tc/Tw  was  1.  A  value  for  Tw  was  assumed  and  the  equations 
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located  in  the  analysis  were  solved  to  obtain  a  secord  value  of  T*.  The  steps 
were  then  repeated  using  thi  s  Tw  as  the  assumed  Tw.  The  process  was  repeated 
until  the  calculated  and  assumed  T^,1  s  agreed  to  within  i*  degrees.  An  upper 
limit  of  30  iterations  was  fed  into  the  program. 

!Ihe  temperature  profile  in  the  laminar  region  was  then  calculated. 

A  biv&rlent  map  for  Eu  ■  1.0  plus  other  similar  maps  for  Eu  *  .5  and 
Eu  w  Overs  then  used  to  obtain  (l  -  Cp  ).  A  value  of  (1  *=*  Q  )  was  obtained 
from  each  of  the  maps  in  the  same  way  a 3  in  the  first  section.  Once  Eu  was 
known,  second  order  interpolation  could  then  be  used  on  the  three  values  to 
obtain  (1  -  )•  T*-  was  iterated  on  just  as  in  the  first  section  except  that 

convergence  is  not  always  obtained  after  30  iterations.  This  occurs  whenever 
fw  is  greater  than  1.  As  this  condition  is  physically  meaningless,  Tw  and 
nr»c  were  printed  as  zero  when  convergence  was  not  attained. 

The  third  segment  of  the  program  deals  with  the  turbulent  flow  region. 
It  utilizes  an  analytic  expression  for  (l  -  (P  )  rather  than  maps.  7^,  is  iter¬ 
ated  on  as  in  the  first  two  steps. 

The  computer  time  for  any  given  set  of  input  parameters  is  about  one 

minute. 


The  sample  data  used  to  "debug"  the  computor  program  are  contained 
in  Table  II. 

The  calculated  results  for  these  data  are  presented  in  Table  III  in 
the  form  that  is  printed  out  by  the  machine.  The  syiribols  printed  by  the  computor 
and  the  corresponding  notation  as  used  in  this  report  are  both  indicated  in  the 
column  headings  of  Table  HI. 
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2  (No«  of  x*s) 

3  (NO. 

of  y's) 

» 

.017  (Xi) 

.155  C:*) 

.0085  (Si) 

/ 

.0035  (I)*) 

350.Cw) 

1000. 

(Tc,r) 

1025. 

GO 

1.  Ov) 

.7(Pr)  1100.  (if)  1' 

5.  C™) 

-10.  (r\) 

(where  K/T-rrv 

•O") 

C 


.017 


.017 


H,*- 

w 

11500. 

0. 

noce. 

io:o. 

10500. 

190c. 

aw 


2J4C0.  ?500.  .0020?  600. 


u5oo. 


.1 55  .0354  11000. 


.1 55  .0356 


19X\ 


23C0, 

10850. 

.00292 

0. 

2J-00, 

9500. 

•  CC209 

60c. 

2300. 

7520. 

.00306 

1780. 

Calculation  identification  number 


Typical  input  data  supplied  to  the  I.B.K.  70U  Computer  in  the  form  used 
by  the  key  puncher. 
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TABULATION  OF  COMPUTED  RESULTS 
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